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Abstract 

A  general  procedure  is  presented  for  computing  the  gas  phase  heats  of  formation  of  a  wide 
y  oret%  of  organic  compounds  A£  for  the  formation  of  the  molecule  from  its  dements  at  0°K  is 
obtained  from  density  functional  calculations  Gaussian  92/DFT)  for  optimized  geometries  This 
resui:  s  .onserted  to  AH  at  by  assuming  ideal  behavior  and  adding  the  translational, 

rotaconai  and  vibrational  energies  Additive  correction  terms  corresponding  to  the  various 
.  s  r  brutn'R  states  of  the  carbons,  nitrogens  and  oxygens  were  developed  using  a  data  base  of  54 
„  'm rounds  The  experimental  AH,-  values  of  these  compounds  are  then  reproduced  with  an 

j  .er.,cc  xbsoiute  error  of  .'  kcaimole  and  a  standard  deviation  of  4  keal/moie  For  a  group  of  ten 
•fs;  .  dH-s  -hat  were  not  part  of  the  data  base,  the  average  absolute  error  is  5  5  kcal/mole  and  the 

standard  de  Mat  ion  is  4  :  kealmoie 
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imcoflufiUon 

In  evalu  ic  potential  usefulness  of  a  proposed  but  not  yet  synthesized  compound,  it 
may  be  of  conquerable  importance  to  have  a  reliable  estimate  of  its  heat  of  formation,  AHf,  since 
this  can  be  a  key  to  predicting  the  compound's  thermochemical  behavior.  A  variety  of  schemes  for 
obtaining  heats  of  formation  theoretically  have  indeed  been  introduced,  based  on  such  concepts  as 
’Mxiesmic  and  homodesmic  reactions,  atom  equivalents,  group  equivalents,  etc.1'12 

A  very  direct  approach,  which  we  use  m  this  work,  is  to  simply  compute  AH  for  the 
tormation  .'t  the  compound  from  its  elements  For  this  purpose,  we  wish  to  exploit  the  advantages 
.if  non-locu:  density  functional  theory  ,  which  permits  energies  to  be  calculated  for  large  systems  at 
a  e%e!  of  accuracy  comparable  to  correlated  ab  initio  methods  ’  *J  Our  objecuve  is  a  general 
procedure  that  can  be  applied  to  a  brr  of  organic  compounds,  and  is  not  limited  to  any 

specific  classes 

Pnxedurc  and  Results 

The  standard  heat  of  formation  of  a  compound.  vHt.  is  define  d  as  AH  for  the  reaction  in 
which  it  is  formed  from  us  elements  in  their  most  stable  forms,  when  they  and  the  compound  are 
all  at  pressures  of  one  bar  Thus.  AH,  for  liquid  nitremethane  at  25s C  refers  to  the  reaction. 

C  (graphite)  +  ^H;(g)  +  ;  N;  (g)  +  CMg)  — ►  CH3NO7  (1)  (1) 

In  our  computational  approach,  however,  we  would  calculate  AH  for  eq.  (2).  convert  it  to  AH.  and 
then  combine  the  latter  with  the  experimental  enthalpy  of  formation  of  gaseous  carbon  and  the 
enthalpy  of  condensation  of  gaseous  nitromethane. 

C(g)  +  |H2(g)  +  j  N2  (g)  +  Ch(g)  ->  CHsNChfg)  (2) 

We  compute  the  optimized  geometries  and  corresponding  energies  and  vibration 
frequencies  for  the  reactants  and  product  in  each  formation  reaction  with  the  density  functional 
option  of  Gaussian  92/DFT,25  using  the  Becke  exchange26  and  the  Perdew-Wang  correlation  27 
functionals  and  a  6-3  lG(d,p)  basis  set.  From  these  energies  is  determined  AE. 

This  AE  corresponds  to  0°K  and  fixed  nuclei.  It  is  converted  to  AE  at  T°K  by  adding  the 
translational,  rotational  and  vibrational  energies,  assuming  ideal  behavior:5 

AEtrans  =  ^RT  (3) 


AE 


rot  - 


|  RT  (non-linear  molecule) 


(4a) 


2 


AErot  =  RT  (linear  molecule)  (4b) 

^«b  =  ihIv.+»»IvJe*pM-l]  ‘  (5) 

i  i  L  V  / 

where  N  is  Avogadro’s  number  and  the  other  constants  have  their  usual  meanings,  v;  is  the 
vibration  frequency  of  the  ith  normal  mode.  Then, 

AH(T)  =  AE(T)  +  ( 1  -nr)RT  (6) 

where  nr  is  the  total  number  of  moles  of  gaseous  reactants  in  the  formation  reaction. 

When  this  approach  was  used  to  compute  the  gas  phase  heats  of  formation  at  25°C  for  a 
group  of  54  molecules  (primarily  organic),  the  results  designated  in  Table  I  as  " uncorrected’’  were 
obtained.  While  some  of  these  compare  well  with  the  experimental  values,  others  are  very  poor. 
The  average  absolute  error  is  22  kcal/mole. 

In  order  to  improve  the  overall  accuracy  of  this  method,  we  introduced  a  set  of  additive 
correction  terms.  Xj  After  experimenting  with  various  options,  we  settled  upon  having  a  Xj  for 
each  coordination  state  of  carbon,  nitrogen  and  oxygen.  Their  values  were  determined  by 
requiring  that  the  calculated  heats  of  formation  reproduce  as  well  as  possible  the  experimental  ones 
in  our  data  base.  Thus  we  minimized  the  quantity  Z. 


where 

ci=Sk«iXJ-Aj  (8) 

i=i 

In  eq.  (8),  the  summation  is  over  the  n  correction  terms,  each  of  which  is  applied  kjj  times  and 
only  on  the  right  side  of  the  formation  reaction  for  molecule  i.  A*  is  the  error  in  its  calculated 

Aj  =  AHf  .(calc)  -  AH*  (exp)  (9) 

gaseous  heat  of  formation,  and  m  is  the  number  of  molecules  in  the  data  base,  54.  The 
minimization  of  Z  was  carried  out  by  means  of  a  multidimensional  pattern  searcn  in  which  all  Xi  arc 
varied  simultaneously  through  ^mall  steps,  the  lengths  of  which  are  adjusted  individually  so  as  to 
provide  the  steepest  descent  of  Z.28-29 


The  resulting  correction  terms  X\  are  listed  in  Table  H  When  these  are  appropriately  added 

O 

to  the  calculated  uncorrected  AHf  values  in  Table  I,  the  "corrected"  AHf  (calc)  are  obtained.  These 
reproduce  the  experimental  values  quite  well;  the  average  absolute  error  and  standard  deviation  are 
3  and  4  kcal/mole,  respectively. 

In  order  to  test  our  procedure,  we  computed  the  gas  phase  heats  of  formation  for  ten 
compounds  that  were  not  part  of  the  data  base  used  to  determine  the  correction  terms,  but  for 
which  experimental  values  are  known.  The  results  are  in  Table  III.  The  average  absolute  error  is 
3.5  kcal/mole  and  the  standard  deviation  is  4. 1  kcal/mole;  these  are  very  similar  to  what  was 
obtained  for  the  compounds  in  the  data  base. 

Conclusion 

We  conclude  that  the  computational  approach  that  has  been  presented  is  an  effective  means 
for  obtaining  reliable  estimates  of  the  gas  phase  heats  of  formation  of  a  wide  variety  of  organic 
compounds.  It  is  anticipated  that  the  scope  and  accuracy  can  be  increased  by  expanding  the  data 
base  used  to  obtain  the  correction  terms. 
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Table  L  Gas  phase  heats  of  formation. _ 

_  Heats  of  formation,  kcal/mole. 


Molecule 

Expt.a 

Calc., 

uncorrected 

Error 

Calc., 

corrected 

Error 

chlorodifluoromethane 

-115.3 

-1 16.14 

-0.8 

-110.46 

4.8 

formic  acid 

-90.51 

-85.91 

4.60 

-84.40 

6.11 

1-fluoropropane 

-68.33 

-82.31 

-13.98 

-65.27 

3.06 

urea 

-58.75 

-61.37 

-2.62 

-50.56 

8.19 

water 

-57.79b 

-42.18 

15.61 

-47.31 

10.48 

ethanol 

-56.21 

-58.77 

-2.56 

-52.54 

3.67 

acetone 

-51.94 

-64.96 

-13.02 

-46.95 

4.99 

methanol 

-48.16 

-44.46 

3.70 

-43.92 

4.24 

dimethyl  ether 

-44.00 

-50.41 

-6.41 

-44.19 

-0.19 

neopentane 

-40.18 

-65.36 

-25.18 

-36.96 

3.22 

n-butane 

-30.02 

-53.17 

-23.15 

-30.45 

-0.43 

methyl  nitrate 

-29.21 

-43.05 

-13.84 

-39.21 

-10.00 

chloroethane 

-26.79 

-39.57 

-12.78 

-28.21 

-1.42 

formaldehyde 

-25.96 

-25.75 

0.21 

-19.10 

6.86 

propane 

-25.02 

-42.35 

-17.33 

-25.31 

-0.29 

nitroethane 

-24.45 

-42.73 

-18.28 

-20.11 

4.34 

ethane 

-20.0 

-31.63 

-11.6 

-20.27 

-0.3 

methane 

-17.8 

-22.56 

-4.8 

-16.88 

0.9 

nitrome  thane 

-17.8 

-28.13 

-10.3 

-11.19 

6.6 

methyl  nitrite 

-15.9 

-27.80 

-11.9 

-19.38 

-3.5 

ammonia 

-11. 0b 

-7.56 

3.4 

-5.48 

5.5 

methylamine 

-5.50 

-11.17 

-5.67 

-3.42 

2.08 

isobutene 

-4.04 

-26.09 

-22.05 

-2.64 

1.40 

oxazole 

-3.70 

-24.26 

-20.56 

-3.99 

-0.29 

dimethy  lmtramxne 

-1.1 

-27.36 

-26.3 

-2.67 

-1.6 

propene 

4.78 

-11.41 

-16.19 

6.37 

1.59 

rutroethylene 

9.c 

-9.14 

-18. 

14.21 

5. 

4,5-dihvdro-3- 

9.32 

-18.84 

-28.16 

11.96 

2.64 

nitroisoxazole 

l,3,5-tria2jne-2,4,6- 

12.4 

-28.24 

-40.6 

17.95 

5.6 

triamine 

chlorobenzene 

12.4 

-26.88 

-39.3 

9.39 

-3.0 

ethylene 

12.5 

3.53 

-9.0 

15.63 

3.1 

cyclopropane 

12.7 

-4.87 

-17.6 

12.17 

-0.5 

1,3-dmitrobenzene 

12.9 

-41.71 

-54.6 

17.09 

4.2 

1 ,4-dinitropiperazine 

13.9 

-37.30 

-51.2 

12.10 

-1.8 

4-mtroaxuline 

14.1 

-36.71 

-50.8 

12.91 

-1.2 

acetonitrile 

15.4 

5.45 

-10.0 

15.37 

0.0 

nitrobenzene 

16.1 

-31.25 

-47.4 

16.29 

0.2 

isoxazole 

18.8 

-4.78 

-23.6 

15.49 

-3.3 

benzene 

19.7 

-17.32 

-37.0 

18.95 

-0.8 

aniline 

20.8 

-19.71 

-40.5 

18.64 

-2.2 

methylhydrazine 

22. 6b 

11.11 

-11.5 

20.95 

-1.7 

dimethylfuroxan 

24.45 

-19.83 

-44.28 

16.40 

-8.05 

dimethylfurazan 

25.65 

-13.11 

-38.76 

19.75 

-5.90 

pyrrole 

25.88 

-1.99 

-27.87 

24.27 

-1.61 

(continued) 
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Table  I.  Gas  phase  heats  of  formation  (continued). 


Heats  of  formation,  kcal/mole. 

Molecule 

Expt.a 

Calc., 

uncorrected 

Error 

Calc., 

corrected 

Error 

1,3-butadiene 

26.29 

4.26 

-22.03 

28.44 

2.15 

indole 

37.40 

-17.17 

-54.57 

33.28 

-4.12 

propyne 

44.19 

31.32 

-12.87 

42.64 

-1.55 

allene 

45.53 

26.95 

-18.58 

45.08 

-0.45 

l,3,5-trinitro-l,3,5- 

45.8d 

-16.96 

-62.8 

40.09 

-5.7 

triazacylohexane 

pyrimidine 

46.82 

8.39 

-38.43 

47.11 

0.29 

pyrazine 

46.87 

12.00 

-34.87 

50.73 

3.86 

acetylene 

54.54 

50.44 

-4.10 

56.09 

1.55 

pyridazine 

66.52 

28.78 

-37.74 

67.51 

0.99 

tetrazoie 

79.95 

49,52 

-30.43 

79.46 

-0.49 

aUnless  otherwise  indicated,  experimental  values  are  from  J.  B.  Pedley,  R.  D.  Naylor  and  S.  P. 
Kirby,  Thermochemical  Data  of  Organic  Compounds ,  2nd  ed..  Chapman  and  Hall:  London,  1986. 
b Handbook  of  Chemistry  and  Physics,  71st  ed„  D.  R.  Lide,  ed.,  CRC  Press:  Boca  Raton,  FL, 
1990.  This  is  also  the  source  for  the  enthalpy  of  formation  of  gaseous  carbon. 

CY.-R.  Luo  and  J.  L.  Holmes,  J.  Phys.  Chem.,  26,  9568  (1992). 
dD.  R.  Stull,  E.  F.  Westrum  and  G.  C.  Sinke,  The  Chemical  Thermodynamics  of  Organic 
Compounds,  Wiley:  New  York,  1969;  J.  M.  Rosen  and  C.  Dickinson,  J.  Chem.  Eng.  Data,  14, 
120  (1969). 
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Table  EL  Calculated  correction  terms. 


Atom 


Coordination  number  Correction  term,  kcal/mole 


C  4  5.68 

C  3  6.05 

C  2  2.82 

N  3  (planar)  10.04 

N  3  (pyramidal)  2.08 

N  2  7.27 

N  1  1.42 

O  2  -5.14 

O  1  0.61 


Table  III.  Gas  phase  heats  of  formation:  test  cases. 


Heats  of  formation,  kcai/moie 

Molecule 

Expt.a 

Calc., 

corrected 

Error 

methyl  acetate 

-98.45 

-97.04 

1.41 

ethylene  glycol 

-92.61 

-85.62 

6.99 

1,4-dioxane 

-75.48 

-79.61 

-4.13 

diethyl  ether 

-60.25 

-63.41 

-3.16 

acetamide 

-56.96 

-53.72 

3.24 

benzaldehyde 

-8.77 

-9.71 

-0.94 

furan 

-8.34 

-12.96 

-4.62 

vinyl  chloride 

8.91 

3.41 

-5.50 

1 ,3 ,5-trinitrobenzene 

14.9 

18.23 

3.3 

pyridine 

33.56 

31.39 

-2.17 

aJ.  B.  Pedley,  R.  D.  Naylor  and  S.  P.  Kirby,  Thermochemical  Data  of  Organic  Compounds , 
2nd  ed..  Chapman  and  Hall:  London,  1986. 


